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ABSTRACT 

Ultracompact Minihalos (UCMHs) have been proposed as a type of dark matter sub¬ 
structure seeded by large-amplitude primordial perturbations and topological defects. 
UCMHs are expected to survive to the present era, allowing constraints to be placed 
on their cosmic abundance using observations within our own Galaxy. Constraints on 
their number density can be linked to conditions in the early universe that impact 
structure formation, such as increased primordial power on small scales, generic weak 
non-Gaussianity, and the presence of cosmic strings. We use new constraints on the 
abundance of UCMHs from pulsar timing to place generalised limits on the parameters 
of each of these cosmological scenarios. At some scales, the limits are the strongest to 
date, exceeding those from dark matter annihilation. Our new limits have the added 
advantage of being independent of the particle nature of dark matter, as they are 
based only on gravitational effects. 

Key words: dark matter, early Universe, cosmological parameters, cosmology: mis¬ 
cellaneous 


1 INTRODUCTION 

Rare objects seeded by large density contrasts at early times 
are an effective tool for probing the early universe. By deter¬ 
mining the present-day abundance of such objects, we can 
investigate the processes leading to their production. Ex¬ 
amining the formation of different objects allows us to test 
different physical processes, scales, and epochs. This probe 
of the early universe is unique in its approach, allowing for 
constraints to be placed on cosmological parameters at scales 
far smaller than would otherwise be accessible. 

Primordial black holes (PBHs) are an extreme example 
of such objects, and have long been used as a probe of the 
small-scale universe. Should a primordial fluctuation exceed 
a large threshold value at horizon entry (<5 > 0.3), the force of 
gravitation will overcome that of pressure and the region will 
collapse, forming a black hole (Carr & Hawking 1974; Carr 
1975). Constraints on the production of primordial black 
holes have been obtained from a multitude of methods, and 
have been used to weakly constrain curvature perturbations 
(Josan et al. 2009; Carr et al. 2010; Alabidi et al. 2012), non- 
Gaussianity (Young & Byrnes 2013; Shandera et al. 2013; 
Young & Byrnes 2015a,b) and inflation (Bringmann et al. 
2002; Peiris & Easther 2008) over a relatively large range of 
scales. 
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In cases where density fluctuations are larger than 

5 ~ 10“^ but too small to form a PBH, the dark mat¬ 
ter contained in the perturbation is expected to collapse 
so quickly that an ultracompact minihalo (UCMH) would 
form (Berezinsky et al. 2003; Ricotti & Gould 2009; Scott 

6 Sivertsson 2009). UCMHs are distinguished from regular 
dark matter structure by their very early time of collapse, 
around the time of matter-radiation equality or even ear¬ 
lier (Berezinsky et al. 2012). Consequently, UCMHs have 
extremely steep density profiles, and are expected to persist 
to the present day, as they would not be easily tidally dis¬ 
rupted (Berezinsky et al. 2006, 2008; Bringmann et al. 2012). 
It has been shown that limits on the abundance of UCMHs 
can be mapped to corresponding limits on processes that 
are expected to increase their production: increased primor¬ 
dial power on small scales (Josan & Green 2010; Bringmann 
et al. 2012), non-Gaussianity (Shandera et al. 2013), and the 
presence of cosmic strings in the early Universe (Berezinsky 
et al. 2011; Anthonisen et al. 2015). UCMHs have also been 
studied extensively for their promise as sources of dark mat¬ 
ter annihilation or decay (Scott & Sivertsson 2009; Lacki & 
Beacom 2010; Yang et al. 2011a,b,c, 2013a,b,c; Zhang 2011; 
Zheng et al. 2014). 

To date, the strongest limits on the UCMH abundance 
have come from non-detection of dark matter clumps in 
gamma rays by the Fermi Large Area Telescope (LAT), 
relying on the assumption that dark matter can annihilate 
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(Bringmann et al. 2012). In Paper I (Clark et al. 2015), 
we showed that a population of UCMHs also will pro¬ 
duce a detectable effect on the period derivative of pulsars, 
due to their gravitational time delay. By exploiting this ef¬ 
fect, we showed that purely gravitational arguments place a 
strong limit on the fraction of dark matter within the Milky 
Way contained within UCMHs. These limits are significantly 
stronger than the only previous gravitational limits (Zack- 
risson et al. 2013), which were placed by assuming non¬ 
detection of small distortions in the images of macrolensed 
quasar jets. Although they do not cover as broad a mass 
range, for some masses the pulsar limits are even stronger 
than those from gamma-ray searches. 

Here we apply our new pulsar limits on the UCMH 
abundance to produce updated, fully model-independent 
constraints on cosmological scenarios that could give rise 
to UCMHs. By calculating the expected UCMH abundance 
for a given primordial power and scale, in Section 2 we give 
generalised constraints on the small-scale primordial power, 
as well as on simple power-law spectra. As the production of 
rare objects has been seen to be very sensitive to higher mo¬ 
ments of the distribution of primordial fluctuations (Bullock 
& Primack 1997; LoVerde et al. 2008; Shandera et al. 2013), 
in Section 3 we place limits on the amount of generic non- 
Gaussianity allowed on small scales. Cosmic strings - topo¬ 
logical defects from symmetry-breaking phase transitions in 
the early Universe (see e.g. Brandenberger 1994) - have also 
been shown to act as seeds for formation of dark matter sub¬ 
structure (Berezinsky et al. 2011; Anthonisen et al. 2015). 
In Section 4, we apply our new limits on the UCMH abun¬ 
dance to constrain the cosmic string tension. In what fol¬ 
lows, unless stated otherwise, we closely follow the methods 
of Bringmann et al. (2012), Shandera et al. (2013) and An¬ 
thonisen et al. (2015) for the respective cosmological scenar¬ 
ios. The code used by each has been implemented in v5.1.2 
of DarkSUSY (Gondolo et al. 2004), providing routines to 
compute the abundance of UCMHs independent of the as¬ 
sumed model of dark matter. 


2 CONSTRAINTS ON PRIMORDIAL POWER 

Primordial fluctuations are thought to have given rise to the 
large-scale structure of the universe. These density pertur¬ 
bations acted as the seeds for small-scale structures, which 
gravitationally collapsed and merged to form a network of 
sheets, filaments, and voids (for an overview see Mo et al. 
2010). These fluctuations are very well constrained on large 
scales by many observations (McDonald et al. 2006; Reid 
et al. 2010; Chluba et al. 2012; Hinshaw et al. 2013; Planck 
Collaboration et al. 2014; Planck Collaboration 2015), indi¬ 
cating that their power appears to be nearly equal on all such 
larger scales. However, given that the imprint of fluctuations 
has yet to actually be observed on small scales, mechanisms 
that would increase (or decrease) power on small scales by 
some degree are not disallowed. 

While the PBH abundance has been used as a probe 
of curvature perturbations on very small scales (as small 
as fc ~ 10^® Mpc“^), this constraint is quite weak in com¬ 
parison to others - by a factor of approximately 7 orders 
of magnitude. Similarly, the UCMH abundance has been 
used to constrain curvature perturbations far more strongly 


{Vti < 10“^ in the range 10 < fc < 10^ Mpc“^), using 
gamma-ray searches with the Fermi-LAT (Bringmann et al. 
2012). From the limits on the present-day UCMH number 
density within the Milky Way that we found in Paper I, here 
we go beyond the assumption of annihilating dark matter 
and provide concrete, model-independent limits on primor¬ 
dial power at small scales. 

The present day mass of a UCMH is related to the 
co-moving radius R of the initial overdense region at horizon 
entry by (Bringmann et al. 2012) 
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The fraction of dark matter expected to be contained 
in UCMHs of mass M5 at redshift z is defined as 




( 2 ) 


where /3(i?) is the probability that such a region will seed the 
formation of a UCMH, z^q is the redshift of matter-radiation 
equality, and is the fraction of matter that 

is cold dark matter (CDM). Here accretion of dark matter 
from the cosmological background onto UCMHs is taken to 
continue up to « ~ 10, after which structure formation has 
evolved such that the majority of halos will be within grav¬ 
itationally bound systems. The present day fraction in the 
Milky Way will then be fiviw ~ 250.77/3(77). 

Assuming that the primordial perturbations follow a 
Gaussian distribution, the probability of UCMH formation 
may be found as 


m) 


o'x,h(-R) 

^ exp 




( 3 ) 


where cr® h is the dark matter mass variance at horizon en¬ 
try and 5^ is the minimum density contrast required to pro¬ 
duce a UCMH. The minimum density contrast is a function 
of both wavenumber k, and the latest allowed redshift of 
UCMH collapse, z^ (see Appendix A in Bringmann et al. 
2012). At redshifts 2 ; < Zc, the background dark matter will 
collapse with sufficient angular momentum that the radial 
infall approximation required to produce the steep density 
profile that characterises a UCMH no longer applies. As 
there is not yet a concrete understanding of what this latest 
redshift of collapse is, we display results for both a rather 
conservative estimate of Zc = 1000 and the slightly more 
liberal Zc = 200. 

Solving Eq. 3 by use of Brent’s Method (Brent 1973), 
in conjunction with the limits on /mw in Paper I, we find 
constraints on primordial mass variance, u® h- To express 
the amplitude of a curvature perturbation Viz in terms of the 
mass variance, a power spectrum model must be assumed. 
We follow the power spectrum normalisation described in 
Appendix B of Bringmann et al. (2012), for three different 
models: 


(i) A ‘generalised’ power spectrum, which assumes local 
scale invariance rather than the global invariance of the 
Harrison-Zel’Dovich model: 

rn{k)=rn{kR)i^—j . (4) 

Here naikii) is the local slope of the power law at kR, which 
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Figure 1. Upper limits (at 95% CL) on the amplitude of pri¬ 
mordial curvature for a generalised power spectrum. We show 
those obtained from both gamma-ray searches and pulsar tim¬ 
ing, and for two assumed latest allowed redshifts of UCMH col¬ 
lapse, 2c = 200 and 2c = 1000. The constraints obtained by 
gamma-ray searches are shown for an assumed dark matter mass 
of = 1 TeV, which annihilates entirely into bb pairs with 
thermally-averaged cross section {(Jv) = 3 x 10“^® crn^s”^. 


Figure 2. Limits on the spectral index, for a scale-free pri¬ 
mordial power spectrum, considering only constraints on ^ 
at wavenumbers smaller than k. These constraints are derived 
from 95% CL upper limits on UCMH number density from both 
pulsar timing and gamma-ray searches, for two redshifts of lat¬ 
est collapse, 2 c. Gamma-ray search limits assume the same dark 
matter model as those in Fig. 1. 


we take to be tir = 1. It should be noted that the limits we 
derive are expected to change for ur ^ 1, as it is not possible 
to relate mass variance and curvature entirely without model 
assumptions. This generalised power spectrum provides a 
normalisation of 

o-x,h(^) = 0.907p7j(fc), (5) 

resulting in limits on primordial curvature Viz. 

We show the resulting limits in Fig. 1. For the case 
Zc. = 200, these limits are of comparable strength to those 
obtained from large-scale observations (log^^g'PTj < —8.5), 
but are extended to much larger k. We reiterate, however, 
that it it is not currently known if UCMH formation can 
continue up to this point, so the weaker limits {zc = 1000) 
should be considered more robust. 

(ii) A scale-free spectrum with constant spectral index Us'. 

'P(fc) oc (6) 

We again follow the method in Bringmann et al. (2012), 
with our derived constraints on spectral index shown in Fig. 
2. The appropriate limit to take from these constraints will 
be the lowest at any scale: Us < 1.24 {zc = 1000), and Us < 
1.02 {zc = 200). Although neither of these constraints is as 
strong as the corresponding limit from gamma-ray searches 
(ris < 1.16 for Zc = 1000, and Us < 1.00 for Zc = 200), they 
apply without any assumptions about the specific particle 
nature of dark matter. Likewise, the limits on the scale-free 
spectral index from cosmic microwave background (CMB) 
observations are in agreement with those we find here (e.g. 
Us = 0.968 ± 0.006; Planck Collaboration 2015), but are 
markedly stronger. 


(iii) A stepped spectrum - scale-free with spectral index 
Us, with the exception of a discontinuous increase in power 
by p at wavenumber ks : 


V{k) oc A:"® ^ 



for k < ks 
for k > ks 


( 7 ) 


In this case, we assume a constant spectral index of 
Us = 0.968 from CMB observations by Planck, including 
la variations allowable by their measurements. We then de¬ 
rive upper limits on the size of the step p as a function of its 
position ks, which are shown in Fig. 3. We find that for steps 
in the region 10°'® < fcs < 10° Mpc“°, the step size must be 
less than a factor of approximately 11 to 18 (zc = 1000) or 
1.5 to 2.6 {zc = 200), depending upon the location of the 
step and the redshift of latest collapse. In contrast to these, 
limits from gamma-ray searches are mostly independent of 
the wavenumber of the step: Pmax < 10 {zc = 1000) and 
Pmax < 1.7 {zc = 200). 

Even if one assumes the most pessimistic case {zc = 1000, 
Us — 0.974, non-annihilating dark matter), the step size 
must be less than a factor of 18.4 at scales larger than 
fc « 2 X 10° Mpc“^. Although the true upper limit is depen¬ 
dent upon both the redshift of latest collapse and the true 
value of Us, our analysis has provided a strong constraint on 
the size of a step in primordial power at far smaller scales 
than previously available, independent of dark matter anni¬ 
hilation. 
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Figure 3. Constraints on the step size for a stepped primordial 
spectrum of primordial fluctuations, as a function of step position. 
Limits correspond to an assumed spectral index of rig = 0.968, as 
well as 1(T deviations allowed by observations from Planck (Planck 
Collaboration 2015). We show limits from both pulsar timing and 
gamma-ray searches, for redshifts of latest collapse Zc = 1000 and 

Zc = 200 . 


3 CONSTRAINTS ON NON-GAUSSIANITY 

Observations of the CMB suggest that the amplitudes of 
the primordial fluctuations follow a Gaussian distribution. 
However, these observations do not possess the sensitivity 
to rule out a distribution that is only approximately Gaus¬ 
sian. Detection of slight departures from Gaussianity would 
provide considerable insight into the nature of the primor¬ 
dial inhomogeneities. If small deviations from Gaussianity 
are present, the probability of larger amplitude primordial 
fluctuations occurring can be increased, acting to boost pro¬ 
duction of rare objects. Number counts of these objects 
have been shown to be sensitive to any level of deviation 
from a Gaussian distribution - the rarer the object, the 
more sensitive it is as a probe of higher moments of non- 
Gaussianity (NG). Gonstraints on the abundances of both 
PBHs and UGMHs have previously been used to constrain 
non-Gaussianity on small-scales (Bullock & Primack 1997; 
Young & Byrnes 2013; Shandera et al. 2013). However, these 
limits are either very weak (PBHs), or depend on the anni¬ 
hilation of dark matter (UGMHs). From constraints on the 
present-day number density of UGMHs in Paper I, we place 
limits on the level of non-Gaussianity at smaller scales than 
accessible via traditional methods. 

Following the method outlined in Shandera et al. (2013), 
we express the level of non-Gaussianity in terms of a depar¬ 
ture from the locally scale-invariant generalised Gaussian 
spectrum discussed in Section 2. To do this, we use the 
model-independent dimensionless skewness, Ms.. This is a 
generalised form of non-Gaussianity, and may be applied to 
any given model. In this manner, most models that give rise 
to non-Gaussian interactions produce a distribution which 


may be expressed as an Edgeworth expansion 


P{u)dv = 


-1/^/2 


s 

n 


1 + ^ ^ Hg + 2rO 


( 8 ) 


(9) 


1 / Mm + 2 

fcm! V(”^ + 2)!, 

where Hniy) are the Hermite polynomials 

rr / \ f 1 \7i ^ —uj’2. 

HnC) = (-1) e —e 

and Mn are the dimensionless moments of the density con¬ 
trast. Here v = ^x/^ x,h{R) is the ‘rareness’ of a fluctuation 
in the limit of a Gaussian spectrum of perturbations. Re¬ 
membering that the mass is proportional to R by Eq. 1, 
this defines r'min = /(j^^h{R) as the minimum rarity re¬ 

quired to seed the formation of a UCMH of a given mass, in 
the Gaussian limit. The second sum in Eq. 8 is over all sets 
of integers {km} (not members of a single set) that satisfy 
the equation 

s = k\ -\- 2/^2 T ''' “t“ skg. (19) 

Each viable set {km} implies a single value of r, defined as 
r = ki k 2 km- (11) 

Here, higher order moments may each be expressed in 
terms of the third moment Ms- We explore two types of 
higher-moment scaling: hierarchical and feeder, each moti¬ 
vated by particle physics (Barnaby & Shandera 2012). Hi¬ 
erarchical scaling results if the non-Gaussianity is generated 
by a single source, such as inflaton self-interactions or cur- 
vaton models. Otherwise, if non-Gaussian fields are coupled 
to the source of the curvature perturbations, then either the 
feeder scaling or a mixed scaling results. Eor the hierarchical 
scaling, this is expressed as 

n — 2 


M„ = n\ 2 ’ 


.-3 fm 

6 


and for the feeder scaling, as 


= 2 " 


Mi 


1/3 


( 12 ) 


(13) 


In terms of Ms, the probability that a fluctuation of 
comoving radius R at time of horizon entry will produce a 
UCMH is then 


= erfc 


Cl-i- 
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s/Ck 
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for the hierarchical scaling, and 
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Table 1. The reference power, logj^Q^^, used for each limit 
on non-Gaussianity from pulsar timing (Fig. 1) and gamma-ray 
searches (Bringmann et al. 2012), for a range of scales, k, and 
redshifts of latest collapse, Zc- 


k (Mpc 

-2c 

logio Rk 
Pulsar Timing 

logic Rk 

Gamma Rays 

1 X IQi 

200 

-7.12 

-7.62 

1 X lO'* 

200 

-8.50 

-8.40 

2 X 10® 

200 

-8.09 

-8.50 

1 X IQi 

1000 

-5.41 

-5.92 

1 X lO'* 

1000 

-6.79 

-6.71 

2 X 10® 

1000 

-6.39 

-6.87 


for the feeder scaling. Here the integers {km}h are the non¬ 
negative solutions to Eq. 10, and {km}f are non-negative 
integers that obey 

s -f 2 = Ski -h 4/t2 -t- • • • “t“ (s “h 2)/i:s. (lb) 

Again following Shandera et al. (2013), we estimate the 
abundance of UCMHs by eqns. 14 and 15. In order to do 
this computationally, we must truncate the series at some 
finite moment of the distribution. We discard terms with 
powers of Ada greater than 16 for the hierarchical scaling, 
corresponding to all terms with s > 17 in Eq. 14, and terms 
with powers of Ads greater than 17 for the feeder scaling, 
corresponding to all terms with s -|- 2r > 18 in Eq. 15. To 
compensate for this level of truncation, we exclude limits 
for which the estimated error can exceed 20 % (for an in- 
depth description of this error analysis, see section 2.3 in 
Shandera et al. 2013). In order to increase the production 
of UCMHs, the upper tail of the distribution must be larger 
compared to the exactly Gaussian case. This will occur for 
any positive value of Ads. Negative values can increase or 
decrease the contribution of the tail, depending upon the 
relative importance of odd and even km , leading to a strong 
dependence on the order at which the series is truncated; for 
this reason we show only limits on positive Ads. 

From the limits on UCMH abundance from both Bring¬ 
mann et al. (2012) and Clark et al. (2015), we place limits on 
non-Gaussianity as shown in Fig. 4 for both hierarchical and 
feeder scaling. We show the limits as a function of the devi¬ 
ation of the Gaussian power Vnik) from the current upper 
limit 'P-^(fe) that we found previously (i.e. Fig. 1). We give 
our adopted reference powers for each wavenumber, UCMH 
search method, and redshift of latest collapse in Table 1. 

On scales of 1 x 10^ and 2 x 10® Mpc“^ we find that the 
limits from Fermi are both marginally stronger and produce 
limits for a wider range of deviations below their Gaussian 
limit, compared to those from pulsar timing. Conversely, we 
find that on the scale of fc = 1 x 10 "^Mpc“^, pulsars provide a 
stronger limit with comparable extent to those from Fermi. 
These differences between the strength and breadth of each 
of these limits are minimal - however, it must again be noted 
that our limits do not rely on the assumption of annihilating 
dark matter, and so rigorously apply to any dark matter 
model equally, modulo considerations of kinetic decoupling 
and its ability to wash out small-scale structure (see e.g. 
Bringmann 2009). 


4 CONSTRAINTS ON COSMIC STRING 
TENSION 


Cosmic strings are topological defects that may have been 
produced in the early universe, present in many models that 
predict symmetry-breaking phase transitions. Their energy 
is confined within long, thin tubes, forming a vast network 
of infinite-length strings - expected to stretch across the ob¬ 
servable universe. When these strings cross one another (or 
indeed themselves), a section can detach, forming a loop. 
The loops oscillate, radiating gravitational waves, and so 
the cosmic string loops lose energy, eventually decaying away 
completely. These loops can gravitationally accrete matter, 
and thus have been shown to act as seeds for UCMH growth 
(Berezinsky et al. 2011). They have a complicated accretion 
history, dependent upon their time of formation and decay. 
Despite this, the number density of UCMHs expected to be 
formed has been predicted for given string loop radius R and 
tension G^i (Anthonisen et al. 2015). Following their method, 
in combination with the constraints on UCMH number den¬ 
sity from Paper I, we compute constraints on cosmic string 
tension as a function of loop radius. 

The number density of UCMHs of a given mass pro¬ 
duced by cosmic strings is strongly dependent upon the evo¬ 
lution of each string. As such, we must treat different evo¬ 
lution scenarios on a case by case basis, in terms of 4 crit¬ 
ical times: time of loop formation {xi), time of loop decay 
(xd), time of latest allowed UCMH collapse (xc), and the 
time at which UCMH accretion ceases (xto), where time is 
parametrised as 




(17) 


With this parametrisation, the redshift of matter- 
radiation equality, Zsq, corresponds to a: = 1. In what fol¬ 
lows, we take the assumption that xto refers to the time after 
which structure formation has progressed sufhciently to al¬ 
low the majority of UCMHs to be within bound structures, 
preventing further accretion from the smooth cosmological 
background (as discussed in Section 2): zto ~ 10, Xto ~ 284. 
Similar to the previous sections, we examine the case of red- 
shifts of latest collapse of both Zc = 1000 and Zc = 200 , 
corresponding to Xc = 3.12 and Xc = 15.54, respectively. 

We follow a one-scale loop model (Vilenkin 1981; Kib¬ 
ble 1985), which describes loops of a given radius as being 
produced together at the same time. The other critical times 
Xi and Xd are then dependent upon the properties of cosmic 
strings formed at each epoch, as: 


Xi = 


Xd = 



(PR 

\jGu teq 


a 


(18) 

(19) 


where a = 1/2 when x < 1, and a = 2/3 when x > 1, 
and a = 0.05, P = 2n, and 7 = IOtt are constants de¬ 
termined from simulations (Vachaspati & Vilenkin 1985; 
Blanco-Pillado et al. 2011). 

The fractional density of UCMHs of a given mass is 
related to the properties of cosmic strings formed at a par¬ 
ticular time by 


dfMW 

iLMl 




( 20 ) 
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Figure 4. Bounds on non-Gaussianity expressed as dimensionless skewness, Ads, as a function of relative Gaussian power, 'P'R.{k)/'P'R.{k)* 
- given at scales near to the extrema of the limits on Gaussian power from pulsar timing in Fig. 1. Dark shaded regions refer to those 
excluded at 95% GL by either gamma-ray searches (grey) or pulsar timing (blue). Light shaded regions correspond to those that would 
be ruled out had the error due to the truncation of the series in eqns. 14 and 15 not been accounted for. Reference power for a range 
of scales, k, and redshifts of latest allowable collapse, Zc, may be found in table 1. 
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where pdm is the present day density of dark matter in 

the universe, G is a constant dependent upon 
formation and decay of the loop. 

the time of 

^ 1 \ 2 + SXd 

C(xi < l,xa < 1) = , 

3 -1- Sxd 

(21) 

C(xi < l,Xd > 1) = 1, 

(22) 


(23) 

^ 6xto-15xi 

C(x.<l,xd<l) 

(24) 


and n{R, t) is the number density (ignoring decay) of loops of 
a given radius at some time t. During radiation domination 
this is 


n{R, t) = Na‘^/3~^t-^R-^, (25) 


and during matter domination 

n{R, t) = (26) 


where N is another constant to be estimated by simulation. 
We take N = 40, as seen by Blanco-Pillado et al. (2011). 

While these calculations assume that each loop is sta¬ 
tionary, cosmic string loops are expected to be formed pos¬ 
sessing relativistic velocities (Blanco-Pillado et al. 2011). 
This non-zero velocity will decrease the efficiency of UCMH 
accretion, as the infall of matter no longer occurs in a 
spherically symmetric manner. This can be simplistically ac¬ 
counted for by assuming that, if a loop were to travel further 
than some distance KR before decaying, a UCMH will not 
be formed. The differential UCMH fraction will then be sup¬ 
pressed by a factor of 5, resulting in: 

d/MW _ ^I^T^GCNo? ^1/2 

dM° SRp^Di^ ’ ^ ’ 

where X — 1 for loops formed after matter-radiation equal¬ 
ity {xi > 1) or X = ateq/{/3R) for those formed before 
{xi < 1), G is the gravitational constant, and 


_ yj-2 j_2 ^Q'^GpDM{ieq)teq 

K, = J^eq^eq ^ 


3/x 


(28) 


For a particular initial velocity Vi of a cosmic string 
loop, the suppression factor will be: 


5 = 


37rl/2(^2)3/2' 


(29) 


where we again follow Anthonisen et al. (2015) by taking 
^.^ 2 ^ 1/2 _ Q 2 ^ which assumes that the loop velocity dis¬ 
tribution follows that of the long strings. This suppression 
factor may be expressed in terms of K by: 


Vi < K 


{a’yGp.y^'^ 


(30) 


By combining Eq. 27 with the constraints on the UCMH 
abundance from gamma-rays and pulsar timing, we are able 
to place limits on the cosmic string tension for a given loop 
radius and suppression factor (parametrised in terms of K). 
We plot these constraints in Fig. 5 for several values of the 
constant K, adopting a redshift of latest collapse of Zc = 
1000. As this derivation depends on the cosmic string scaling 
solution, it is important to note that the actual resulting 
limit on G/i corresponds to the strongest limit at any R. As 



R (kpc) 


Figure 5. Constaints on cosmic string tension Gfj, as a func¬ 
tion of loop radius, for a range of different velocity suppression 
factors K. We display limits derived from those on UCMH num¬ 
ber density from both pulsar timing (solid lines) and gamma-ray 
searches (dashed lines). Limits from gamma-ray searches again 
assume dark matter mass = 1 TeV, and 100% annihilation 
into bb pairs with cross-section {av) = 3 X 10“^®cm^.s“'^. 


such, not only do the constraints due to pulsar timing remove 
the assumption that dark matter must annihilate, but also 
strengthen the overall best constraint on string tension from 
Gp, < 3.14 X 10“^ to G/i < 6.49 x 10“® for loops that are 
able to travel 1000 times their own radius and still form a 
UCMH. 

The shape of these limits may be understood by com¬ 
parison with the constraints on UCMH number density. For 
subhalo masses greater than approximately W^Mq, these 
limits are saturated by the probability that there are insuf- 
hcient UCMHs within the Milky Way to provide a reliable 
signal. The strongest limit at this point is mapped to a con¬ 
straint on the mass of the loop, which in turn is proportional 
to its tension and radius as M = pj3R. This mapping trans¬ 
lates the peaked shape of constraints on UCMH abundance 
directly to those found here. 


5 CONCLUSION 

The large-scale structure of the universe is thought to have 
been seeded by small fluctuations in the early universe. Al¬ 
though these structures were formed from overdensities of 
order 5 ~ 10~®, larger amplitude fluctuations are thought 
to be able to produce rare structures such as primordial 
black holes (<5 > 0.3) and ultracompact minihalos (0.3 > 

5 > 10 “®). 

The abundance of these rare objects has been used to 
constrain a variety of processes that would boost their pro¬ 
duction beyond that expected from the standard Harrison- 
Zel’dovich (scale-free) model. To date, three such processes 
have been constrained: an increase in primordial power at 
small scales, deviations of the distribution of primordial fluc- 
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tuations from a Gaussian, and the presence of cosmic strings 
in the early universe. By consideration of each process in 
turn, it is possible to link their properties to the present day 
number density of each rare object. 

While previous studies have given constraints on both 
UCMH and PBH number densities, these are either very 
weak, or strongly dependent on assumptions about the spe¬ 
cific particle nature of dark matter. By considering new 
upper limits on the number density of UCMHs (Paper I), 
we provide updated constraints on the properties of each of 
these processes. Although the limits are strongly dependent 
upon the assumed redshift of latest formation of a UCMH, 
Zc, even a very conservative value of Zc = 1000 results in 
some of the strongest limits to date. 

Here we have calculated the contribution of increased 
primordial power on small scales for 3 different power spec¬ 
trum models. For a ‘generalised’ power spectrum, we found 
log^gP-Tj < —6.5 (zc = 1000) and loglo'F’7^ < —8.5 (zc = 
200), in the range 10^ ^ ^ 10^- This is comparable 

to the limits from non-detection of dark matter sources 
by Fermi-LAT. We additionally find limits of Us < 1.24 
{zc = 1000) and ris < 1.02 {zc = 200) on the spectral index 
of a scale-free power spectrum. For a stepped spectrum, the 
non-observation of UCMHs limits the step size to a factor 
of approximately 11-18 (jZc = 1000) and 1.5-2.6 {zc = 200). 

We also provide limits on the dimensionless skewness 
Ais, dependent upon both the scale of the fluctuation, k, 
and redshift of latest collapse, 2c. Assuming two different 
models of scaling with higher moments, we find limits that 
are independent of dark matter annihilation, and (depend¬ 
ing on the scale) are able to be applied to lower primor¬ 
dial power, and to wider variations. Depending on the na¬ 
ture of the non-Gaussianity, these limits can be to be eas¬ 
ily mapped to the more model-dependent quantity /nl, for 
comparison. For example, if we take simple non-linear cou¬ 
pling, TZ{x) = TZa{x) + |/nl \Rg{x)'^ - (7?.G(a:)^)], the two 
may be related by /nl ~ MzIVlfP ■ This results in a con¬ 
straint of /nl < 0(10^) to 0(10®), depending upon the exact 
shape of the primordial power spectrum. 

Finally, we constrain cosmic string tension as G^i < 
6.49 X 10“®, under the assumption that a loop can move up 
to K = 1000 times its own radius and still form a UCMH. 
Although this constraint is stronger than that from CMB 
observations (G/r < 1.7 x 10“^ from Dvorkin et al. 2011), 
this assumed value of K is probably overly optimistic. As 
K decreases, the limit grows significantly weaker. While far 
stronger limits of G/r < 2.8 x 10“® were obtained by Blanco- 
Pillado et al. (2014), their result relies on the proper un¬ 
derstanding of emission of gravitational waves from cosmic 
string cusps. As this is a poorly understood process, these 
constraints must be treated with appropriate caution. 

We have shown that these limits depend heavily on the 
latest redshift at which UCMHs are assumed to be able to 
form. We have displayed constraints for both the conserva¬ 
tive value of Zc. = 1000 and the more optimistic Zc = 200. It 
is important to note that there is, as yet, no strong evidence 
in favour of either value. Should further research be under¬ 
taken to investigate the physical value of Zc, the limits here 
could potentially be improved substantially, leading to the 
exclusion of multiple cosmological models. 
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